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ABSTRACT 

This paper reviews Lilley’s reformulation of Lighthill’s equation and shows that it can be 
obtained as a special case of a much more general acoustic analogy. It also shows how this 
generalized analogy can be used to eliminate some of the difficulties that arise when more 
conventional parallel flow analogies are applied to high speed jets. And, finally, some recent 
applications of these ideas are discussed. 


1. INTRODUCTION 

It’s with considerable pleasure that I write this paper in honor of Geoffrey Lilley 
who is so unassuming that it is easy to forget how important his contributions to 
Aeroacoustics actually are. After thinking about this for awhile, it occurred to me that 
the best way to acknowledge the importance of those contributions would be amplify 
some of his results by embedding them in a more general context. So I started with his 
clarification of the so called “entropy” source term d-ip-c^p) in the Lighthill stress 
tensor 


v = -p'-v’+syp-^pi+o,, 


(i.i) 


that appears on the right hand side of Lighthill’s (1952, 1954) equation 

d 2 p ' 2 d 2 p' = d 2 

dt 2 ° dxjdxj dxfixj l] 


( 1 . 2 ) 


where t denotes the time, x = {x ; } denote Cartesian coordinates, p denotes the density, 
v. denotes the fluid velocity, p denotes the pressure, G / - denotes the viscous stress, the 
subscript zero refers to the uniform background state at infinity, with c {) being its sound 
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speed and p' = p - p 0 being the density fluctuation. Lilley (1996, 2003) introduced a 
modification of Lighthill’s equation that separates this term into isentropic and non- 
isentropic components, which is important because the non-isentropic component is 
proportional to the heat flux vector q. and the viscous work term <3 v. that are usually 
neglected at the high Reynolds numbers where jet noise becomes important. 

While there are other ways of accomplishing this that lead to simpler results 
(Goldstein, 2001) an important advantage of Lilley’s approach is that it can be extended 
in a natural way to a more general formulation that promises to become increasingly 
important as computational capability continues to increase (Bodony and Lele, 2002, 
2003, 2005, Karabasov et al, 2008). The present paper shows how Lilley’s revision of 
Lighthill’s equation can be obtained as a special case of this more general result. Some 
recent applications of this generalized approach are also discussed. 


2. LILLEY’S MODIFICATION OF LIGHTHILL’S EQUATION 

Lilley began by rewriting Lighthill’s equation in terms of the pressure to obtain 


He then noted that the energy equation can be written in terms of the stagnation 
enthalpy 



( 2 . 1 ) 


f , 1 2 
4 = h + — v 

A O 


2 


( 2 . 2 ) 


as 



(2.3) 


where h denotes the enthalpy. The continuity equation 



(2.4) 


can be used to write this as 



(2.5) 
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where h 0 has the obvious meaning. But since h 0 is equal to Cg/(y - 1), where y is the 
specific heat ratio, and p/z - p is equal to h^p/c^ for an ideal (and perfect) gas this 
equation can be written as 


d_ 

dt 


^-co 2 p+(y-i)p v2 / 2 

a 

_i 7 

(k-k) 


" °ij v j ) 

4 J 

+ 3 P V *' 

OXj 

v K j 

=---\ 

h 0 dx t 


( 2 . 6 ) 


in that case. A flow will be isentropic, i.e., the convective derivative of its entropy will 
vanish (Yih, 1977, p. 255), when the viscous and heat conduction effects (i.e., the 
diffusive effects) are zero. So the right hand side of [2.6] will be equal to zero and 
the left hand side will vanish in that case. Equation [2.6], therefore, implies that 
d[{p-c 2 oP)/cl]/dt can be divided into isentropic and non-isentropic (or in Lilley’s 
words diffusive and non-diffusive) components. And inserting it into [2.1] shows that 
(Lilley, 1996, 2003) 

1 d 2 p d 2 p _ d 2 / 

c\ dt 2 dXjdXj dxfixj ^ 1 ; lJ 


d 2 

pv i {h s -h 0 ) + q i -O ij v j 

d 2 

'(y-l)pv 2 ' 

dtdXj 

K 

dt 2 

1 

CN O 

O 

<N 

i 


The advantage of this equation over the original Lighthill result is that the diffusive 
(2 nd term on right hand side) and non-diffusive (3 rd term on right hand side) 
components of the source term are now explicitly displayed. The main disadvantage is 
that Lighthill’ s simple quadrupole source is now replaced by a combination of 
quadrupole, dipole and monopole sources. Lighthill (1954) emphasized the importance 
of displaying the correct multipole order of the source term. He argued that the 
inclusion of extraneous lower order multipole components that are more efficient than 
the higher order multipole components at low frequencies could lead to an over 
estimation of the radiated sound in the low frequency limit. Inclusion of the dipole 
term appears to be unavoidable, but the monopole term can be eliminated by 
introducing the new dependent variable 


P e = P + 



( 2 . 8 ) 


to obtain the somewhat simpler result 


1 <T-p c d 2 p e | 1 9 2 e 4 . 

cl dt 2 dxjdxj dxfiXj cl dtdxj 


(2.9) 
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where the source strength e h is now given by the 4x3 dimensional “Lilley’s Stress 
Tensor” 




Y-l~ 

+ ^ 5 , 


pV +G X 


for X = 1,2, 3, 4 


( 2 . 10 ) 


with 


V 4 = (y-l)(h s -h 0 ), and a 4( - = -(y-l)(^. -c,.,.v.) (2.11) 

(The Greek indices X and k will always range from 1 to 4, all other Greek indices 
will range from 1 to 5 and Latin indices will range from 1 to 3.). 

Lilley dealt with the monopole term by arguing that the relatively small (y— 1 )/ 2 
multiplicative factor should cause it to be negligibly small. But Goldstein and Leib 
(2008) found that it actually makes a significant contribution to the 90° sound field when 
the turbulence is assumed to be quasi-normal and axysymmetric and this probably 
remains true even when it is not quasi-normal. It is now generally believed that the 
remaining enthalpy term (d 2 e Aj /dtdx\ / c\ is an important sound source in heated jets 
and significant effort is underway to develop appropriate models for this quantity 
(Khavaran and Kenzakowski, 2007). 

Equation [2.9] can be made to resemble Lighthill’s equation even more closely by 
introducing the new dependent variable 

X 4 = —Cgt (2.12) 


to obtain 


1 d 2 P e <?Pe = ( 2 . 13 ) 

c\ dt 2 dXjdXj dx x dxj 

where the summation convention is now being used for the Greek as well as the Latin 
indices. The non-linear dependent variable p e (or equivalently p e - p {} ) causes no particular 
difficulty here because it reduces to the ordinary pressure fluctuation p - p {) in the far field 
where the sound is to be calculated. 

3. A GENERALIZED APPROACH 

Lilley was also one of the first to emphasize the importance of removing all non-sound 
generating terms from the source function in order to obtain accurate predictions of the 
sound field. This is very difficult to do in practice, but can be approached more and more 
closely with increased computational effort. The appropriate acoustic analogy equations 
can be obtained by moving more of these silent terms to the left hand side of the Navier- 
Stokes equations and expressing the result as a set of formally linear equations with non- 
linear dependent variables and source terms. 
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This can be accomplished by dividing the dependent variables in the Navier Stokes 
equations, i.e. the density, pressure, velocity and enthalpy into, say, a “base flow” 
component, which is denoted by the over-bars and tildes and “residual” components, 
which are essentially defined by the equations 

p = p + p', p = p + p\ h = h + h' v. = v. + v', (3.1) 

and requiring that the base flow components satisfy the usual hyperbolic conservation 
laws, which can be written fairly compactly as (Goldstein and Leib, 2008) 

D o p = 0, (3.2) 


D p,7 + £ = Ae 

‘ a.t, dxj 11 




y-l 2 


dPe 
dt dx. 


■i V 


y-l 


4 


) 


(3.3) 


(3.4) 


where the base flow variables, as well as the original Navier Stokes variables are assumed 
to satisfy an ideal gas law equation of state, p e denotes a pressure-like variable that can 
differ from the thermodynamic pressure p , e^. denotes an arbitrary 4x3 dimensional 
stress tensor and 


dj ‘% <3 - 5 » 

for any function /. These equations include, among other things, the Euler equations, 
the Navier Stokes equations themselves, and even the Favre filtered Navier Stokes 
equations with very general space-time filtering (Aldama, 1990, Goldstein, 2005, 2008). 

It can now be shown (Goldstein, 2000, 2002, 2003, see also Goldstein and Leib, 2008) 
that the remaining residual variables are determined by the five formally linear equations 

+ (3.6) 


dv, 

Z> 0 «,.+« y — + 


7 V J -'< 


^P'e+P% 

OX; 


= dx^ 


(3.7) 


and 




dv. 

OX- 


V J 




(3.8) 
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where 


C 2 =y p/p 


is the square of the base flow sound speed and 


X ; = 


a _ a 

— + V ■ — 

ax ; dy. 


-ss-(vv.) 


P 9y, 


(3.9) 


(3.10) 


is the base flow advection vector. 

The Navier-Stokes equations have been put into this linearized hyperbolic 
conservation law form by introducing the non-linear dependent variables 


P'e =P-Pe+^Y-P v ' 2 =(p' + -^-W 2 j + P-P e (3.11) 

Uj = pvf, (3.12) 

as well as in the quadratically non-linear source strengths 

e"^e' Xi -e Xi X = 1,2,3, 4 (3.13) 

where denotes the generalized “Lilley Stress Tensor” 

e xi = "P v I v / ,+ ~ _ ^P y,2 + a ?i/ (3.14) 

with 

< = (t-0 /? v- = -(Y-l)(?,.-Oyv;) (3.15) 


and 


h'=h' + -v' 2 (3.16) 

* 2 

is the residual stagnation enthalpy (i.e. the difference between the stagnation enthalpy 
and the base flow stagnation enthalpy) in a reference frame moving with the local base 
flow velocity v .. 
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The nonlinear variable [3.12] is crucial but, as in the Lilley formalism, the nonlinear 
variable [3.11] is mainly introduced to simplify the result. Since the velocity variable is 
already non-linear there is no reason to preclude this potentially very significant (see 
Goldstein and Leib, 2008) reduction in complexity by artificially requiring that the 
pressure variable be linear. The moving frame residual stagnation enthalpy h' s was 
introduced because order of magnitude estimates (Goldstein and Leib, 2008) suggest 
that v 4 = ( 7 - 1 ) h' s will be negligible compared to v\ in cold (i.e. unheated) jets at 
0(1) Mach numbers — which means that e 4[ will then be negligible compared to e and 
that the source term can, therefore, be considerably simplified. But recent RANS 
computotions (Khavaran and Kenzakowski, 2007) suggest that it may be the fixed 
frame residual stagnation enthalpy h ' s + v .vO that is negligibly small compared to v'. in 
cold (i.e. unheated) jets at very high Mach numbers. It would then be better to define v' 4 
by v' 4 =(Y-l)(/7 + v vj ) and, as a consequence, change the right side of [3.8] to 
de" / dxj — (y — 1 ) ^de" / dx j . 

But at present the only safe course of action would be to retain the 4 th component e" 4i 
even for cold jets — in which case the precise definition of v' 4 would be largely irrelevant. 

Equations [3.6] to [3.8] can be written more compactly as 

L w u v =DZ je " for (i,v = 1,2, ...5 (3.17) 

where L^ v and D are simple first order linear operators (given explicitly in Goldstein, 
2006) and u v denotes the five dimensional solution vector 

kM«.^ P '} (3.i8) 

Interpretation of the right side of these equations as acoustic sources implies that that 
the left sides describe the interaction of this sound with as well as its propagation 
through the base flow. The formal linearity of the equations can now be exploited to 
separate out these base flow interaction effects (including the propagation of the sound 
through that flow) from the unsteady source fluctuations that produce the sound by 
using the 4 th component of the adjoint vector Green’s function g^iy,^ \x,t), which is 
related to the 4 th component of the ordinary vector Green’s function g AV (x,t\y,x) by the 
reciprocity relation g“ 4 (y, T\x,t) = g 4V (x,t\y,%) and satisfies the adjoint equation 

K ) yT g“ 4 (*T|*,f) = S^ 4 8(x-y)8(t - r) (3.19) 

(Morse & Fesbach, 1953, p. 870, Tam & Auriault, 1998) where L^ v is the adjoint of L^ y 
to obtain the following formula for the pressure like variable p' e 

p'( x ,t) = -j j y Xj (y,x\x,t)e" j (y,x)dydx 


(3.20) 
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where the propagator (y,r\x,t) is related to g“ 4 (y,r\x,t) by 




*7 


-(y-1)8 m 


d >', 


(3.21) 


p’ e (x,t) reduces to the ordinary pressure fluctuation p\x,t) when x is in the far field and 
I am assuming, for simplicity, that all solid surface effect can be neglected. (But the 
result remains valid in the presence of fixed solid surfaces when viscous effects are 
neglected.) 

Equation [3.17] is similar to [2.9] in that the entire source strength is described by a 
4x3 dimensional tensor. In fact equations [3.2] to [3.4] reduce to the Euler equations 
(with p e = p) when the base flow source strengths e^. are set equal to zero. And these 
equations posses a family of solutions, which quite conveniently provides a good 
approximation to the actual mean flow field in boundary layers, turbulent air jets and 
other nearly parallel flows: namely the unidirectional transversely sheared mean flows 


v t = d n U(x 2 ,x 3 ), p e = constant, p = p(x 2 ,x 3 ). (3.22) 

where U and p can be any functions of the indicated arguments. Equation [3.6], (i.e. the 
continuity equation) then decouples from the remaining four equations [3.7] and [3.8], 
which now become the inhomogeneous compressible Rayleigh equations (Betchov and 
Criminale 1967) 


where 


D o Ui dU dp' e d , 


D„p' dc 2 u, de' iW 3 U 


Dt dxj dxj 


2 * 4 +l ,± 

Dt dt dx x 


(3.23) 


(3.24) 


(3.25) 


denotes the convective derivative. 

It is now well known (Goldstein, 1976) that the velocity variable u [ can be eliminated 
between these equations to obtain the inhomogeneous equation 




dx ; 


dx u 


-— 2 ? S < 


d: 


dx f dx x dx. Dt 2 


d <j ( \ , dU 

ay +h - 1)e '-s;_ 


(3.26) 
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where 



_3_ 

Ox,. 


Dt 2 J dxj dx x dxj 


(3.27) 


is the variable-density Pridmore-Brown (1958) or Rayleigh (Betchov and Criminale, 
1967) operator. This is a form of Lilley’s famous (1974) equation whose source term, 
as emphasized by Morris (personal communication), turns out to be somewhat more 
complicated than the original Rayleigh equation sources. It reduces exactly to equation 
[2.9] in the Lighthill limit 


U — > 0, c — > Cq 


(3.28) 


(after performing a trivial integration with respect to time). The analogue of [2.13] 
would be 




A 




dx 


j y 


OX ■ 


AfL 

dx { dxj 


y-i 


d: 


c e n 


(3.29) 


where 


A i = X = 1,2,3 

OXj 




c 2 Dt ’ 


A = 4 


(3.30) 


While it can be argued that Lilley’s equation (along with its variants) has served as 
the basis of most of the theoretical jet noise research over the past couple of decades, 
this result does have a number of issues. First of all the mean flow interaction effects 
have not been completely eliminated from the source function because v- still has a 
mean flow component. And secondly the source function contains a steady term that 
does not generate any sound and therefore does not belong in the source term. Another 
more serious difficulty is that the Rayleigh operator can support linear instability waves 
that grow without bound in this parallel mean flow model. But these homogeneous 
solutions contribute to the adjoint Greens’ function and therefore to the propagator 
Y,, . {y,x\x,t) when causality is imposed (Briggs,1964, Bers, 1975, Goldstein and Leib, 
2005) — causing the solution to become infinite far downstream in the flow. 

Arguments against imposing causality (Mani 1976; Dowling, Ffowcs Williams & 
Goldstein 1978) usually amount to asserting that it is unnecessary because it is not 
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Figure 1: Jet flow configuration. 


possible to identify an initial time before which the fluctuations (about the base flow) 
have been ‘switched on’ in an acoustic analogy approach and that a boundedness 
requirement can, therefore, be imposed on the solution. But the real problem that needs 
to be solved is not that of a doubly infinite parallel jet but that of a jet issuing from a 
nozzle as shown in figure 1. The steady state (i.e., time stationary) solution to this 
problem can be made unique without considering initial conditions (i.e., causality) by 
imposing a Kutta condition at the nozzle lip. 

Crighton and Leppington (1974) constructed the two dimensional Green’s function 
for the analogous problem of a vortex sheet behind a semi-infinite plate, separating an 
ambient fluid from a stream moving with uniform speed (see also Crighton, 1985, p. 436). 
They found that the imposition of the Kutta condition implies that causality will be satisfied 
in this case. 

The doubly infinite parallel flow model amounts to neglecting the reflection and 
diffraction contributions to the Green’s function in the true jet noise problem so that the 
remaining portion of the solution is likely to be the causal Greens function for the 
doubly infinite parallel flow. This suggest that the parallel flow Green’s function should 
be required to satisfy causality in the jet noise problem — which means that it will not 
remain bounded. 

Another serious difficulty at supersonic speeds is that the adjoint Greens function 

has a singularity at the so called critical level which produces an even stronger non- 
integrable singularity in the propagator y (y,T \x,t) when the observation point x is 
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in the far field — leading to the absurd conclusion that the radiated sound is 
infinitely loud. 

k. THE MEAN FLOW ANALOGY 

None of the difficulties would occur when if the base flow were taken to be the mean 
or time average of the actual flow — even when that flow is a nearly parallel shear flow. 
Because, for example, its slow divergence would cause the instability wave to eventually 
decay instead of growing without bound. And, more importantly at supersonic speeds, 
the flow would no longer posses a critical layer and would, therefore, not have a critical 
layer singularity. 

This base flow is obtained be interpreting the over-bars in equations [3.2] to [3.4] as 
the time averages 




1 T 


(4.1) 


and the tildes as the Favre averages 


• = (p*)/p 

in which case 

P'e=P' + ^Y~( PV' 2 -PV' 2 ) 


(4.2) 


(4.3) 


and more importantly 


= c 


(4.4) 


So the complete residual equation source strength is simply 


e —e 

V] VJ 


(4.5) 


which means it has zero time average and therefore contains no steady component. The 
base flow equations are of course the usual Reynolds averaged Navier-Stokes equations. 
This form of the analogy completely removes the mean flow interaction effects from the 
source term, which is expected to significantly improve the medium and lower 
frequency predictions, but no steady base flow analogy can provide the requisite 
structure for inputting the long range turbulent scattering effects, which tend to become 
more important at the higher frequencies (and could eventually overwhelm the mean 
flow interaction effects at sufficiently high frequencies). 
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But a major advantage of this analogy is that the Green’s function solution [3.20] can 
be used to express the far field pressure autocovariance 

I i 

P 2 {x,t 0 )=— j p'(x,t)p'(x,t + t 0 )dt, (4.6) 

as the convolution product 

p 2 (x,t) = J JJ y KjU (x\y;T],t+t)R KjXl (y;T],t)dy dr] dr (4.7) 

-oo v 

of the “ convolved propagator” 


yK j xi( x \y’ T ld + r)=jy Kj (x\y,t 1 +t + T)y xl (x\y + T],t 1 )dt 1 , (4.8) 


with the two point time-delayed correlation of the source strength 


1 1 

R KjU e"j (y, T 0 )e" (y + T], t 0 +T )dt 0 


(4.9) 


which is, in turn, related to the generalized Reynolds stress autocovariance tensor 
T 

*\jK i{y^)=-^ J [pviv' -pv[v'](j,T 0 )[pv;v;-pv>;](j+v,t 0 +t)fi?t 0 (410) 

-T 


by the relatively simple linear transformation 

(4.11) 


(4.12) 

when viscous effects are neglected. This result shows that there is an exact relation 
between the quantities that are actually measured in Aeroacoustics experiments: namely 
the far field pressure autocovariance and the generalized Reynolds stress autocovariance 
tensor. It is very general and applies to any localized turbulent flow even in the presence 
of fixed solid surfaces. The hope is that most of the non-local “propagation effects” 


jXl ~ j,o m y n^Xl,yn 


where 


= 6 8 


jm 


Y \ ? ? 

2 fiom 
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which would be very difficult to distinguish from the turbulent fluctuations — and, 
therefore, very difficult to model — have been removed from these stresses. 

5. APPLICATIONS OF THE MEAN FLOW ANALOGY 

Karbasov et al (2008) used this form of the analogy to predict the sound from a cold 
Mach 0.75 jet and obtained reasonably good agreement with experimental data from the 
JEAN project. But the Green’s function and, therefore, the propagator is much more 
expensive to compute (in fact prohibitively so for many purposes) than the 
corresponding Lilley equation propagator. There is, however, a way to obtain the best of 
both worlds, by using a perturbation approach that takes advantage of the small spread 
rate, say £, of the jet or shear layer (Goldstein & Leib, 2005, 2008). The base flow 
velocities will then expand like 


v l =U(Y,y T ) + eU il) (Y,y T ) + ... (5.1) 

v T =eV(Y,y T ) + £ 2 V m (Y,y T ) + ... (5.2) 

with similar expansions for the other dependent variables. 

Y^ey, (5.3) 

denotes a slow streamwise source variable that varies on the streamwise length scale of 
the mean flow (see figure 1), and 


y T ={y2^} 

(5.4) 

II 

(5.5) 


denote cross flow variables. I am assuming, of course, that all lengths have been 
normalized by some characteristic cross flow dimension of the jet and all velocities by 
some appropriate characteristic streamwise velocity with similar obvious normalization 
for the density, pressure and temperature. The lowest order terms in this expansion 
correspond to the unidirectional transversely sheared flow with the slow streamwise 
variable entering only parametrically. 

The 4 th component of the adjoint Green’s function is expected to possess a 
corresponding expansion of the form 


8v 4 = 8vf ( y, Y,r\x,t) + £g“i (y,Y, t | x,t) + . . . 


(5.6) 
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where the lowest order term satisfies the adjoint Rayleigh’s equation in the vicinity of 
the source point. But this solution is non-uniformly valid because it still has a critical 
level singularity and becomes unbounded far downstream in the flow. However, the 
method of multiple scales, the WKBJ method and the method of matched asymptotic 
expansions can be combined to obtain a uniformly valid “composite” solution that 
remains finite everywhere in the flow and is still not much more expensive to compute 
than the original parallel flow solution. The resulting “propagator” consists of two 
components. The first of which corresponds to the usual non-causal Lilley equation 
solution appropriately corrected to eliminate the critical layer singularity. The second 
(causal component) is associated with linear instability waves, but differs from 
conventional instability models in that the instabilities, which are now continuously 
‘generated’ along the length of the jet, do not constitute separate sound sources. 

Previous attempts to incorporate instability waves into jet noise predictions have 
argued that there are two distinct sound sources: one associated with the smaller scale 
random motions and the other associated with the large scale coherent structures — with 
the latter source usually being represented by the linear instability waves. But recent 
studies show that the large scale motions are not well represented by these instabilities. 
The present approach avoids this conundrum because the instabilities only appear in the 
propagator with their sole function being to produce the appropriate cause-effect 
relation between the sound and its turbulent source. They do not correspond to any 
physical structures that actually appear in the flow. 

Figure 2 (which is a compilation of data presented in Bridges et al) shows that the 
OASPL directivity patterns are smooth functions of angle at subsonic acoustic Mach 
numbers but exhibit a relatively abrupt increase in slope at about 60° from the 
downstream axis at supersonic Mach numbers — suggesting that some new effect may 
be coming into play at the smaller angles. 

The non-causal component of the Green’s function is expected to produce a relatively 
smooth directivity pattern, but results presented in Goldstein and Leib (2005) show that 
the causal component of the far field Green’s function is highly damped at subsonic 
speeds and also at supersonic speeds when the observation angle (measured from the 
downstream axis) is sufficiently large. However the damping rapidly drops to zero when 
the observation angle become less than 60° (see figures 4 through 6 of Goldstein and 
Leib, 2005) allowing the causal component of the Green’s function to become much 
larger than the non-causal component. But source non-compactness effects can also 
produce an abrupt increase in shallow angle directivity. And Goldstein and Leib (2008) 
tentatively concluded that the causal component only becomes significant when the 
source point is close to the nozzle where reflection and diffraction effects are 
important and can, therefore, be neglected since the latter effects are already assumed 
to be negligible in the weakly non-parallel flow approach. They then used the 
remaining non-causal component (appropriately corrected for the critical layer 
singularity) to predict the sound from cold high speed jets over a range of acoustic 


OASPL(dB) OASPL(dB) OASPL(dB) 


oeroacoustics volume 9 • number 4 <fr 5 • 2010 


415 



40 60 80 100 120 140 160 180 

Inlet angle (deg) 



40 60 80 100 120 140 160 180 

Inlet angle (deg) 



40 60 80 100 120 140 160 180 

Inlet angle (deg) 



40 60 80 100 120 140 160 180 

Inlet angle (deg) 



40 60 80 100 120 140 160 180 

Inlet angle (deg) 



40 60 80 100 120 140 160 180 

Inlet angle (deg) 


Figure 2: OASPL directivity patterns measured by Bridges on the Glenn SHJAR 

rig (Bridges et al, 2009). 
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Mach numbers from 0.5 to 1.4 and were able to obtain reasonable agreement with 
measurements taken on the Glenn SHAR rig (Khavaran, Bridges & Freund, 2002). 
These issues obviously need additional consideration 

6. CONCLUSIONS 

Lilley developed a modification of Lighthill’s equation that explicitly separates the so 
called entropy source into its isentropic and non-isentropic components. The present 
paper shows how this result can be obtained as a special case of a generalized acoustic 
analogy that can be used to eliminate some of the difficulties that arise when more 
conventional parallel flow analogies are applied to high speed jets. 
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